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Abstract 

Necessary  and  sufficient  conditions  for  the 
exponential  stability  of  equilibria  are  obtained  for 
the  resistive^  viscous,  incompressible  MHD  sheet  pinch 
in  a  gravitational  field.   Insulating  as  well  as  conduct- 
ing boundaries  are  considered.   The  unperturbed  magnetic 
field  B^  (horizontal),  resistivity  t)  ,  viscosity  v„, 
and  mass  density  p„  are  pennitted  to  be  arbitrary 
(nonconstant )  functions  of  the  vertical  coordinate 
consistent  with  the  equilibrium  equations.   In  particular, 
it  is  shown  that  if  rj   is  not  constant  and  the  boijndaries 
are  insulating,  the  pinch  cannot  be  gravitationally 
stabilized. 
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1 .   Introduction. 

Resistive  plasma  instabilities  have  received 
considerable  study  in  recent  years.   In  particular, 

instabilities  in  the  resistive  incompressible  MHD  sheet 

1  -7 
pinch  have  been  investigated  by  several  authors.   ' 

Most  of  this  work  has  been  restricted  to  the  limiting  cases 

of  large  or  small  resistivity."^'    Roberts  and  Taylor 

considered  purely  gravitational  resistive  instabilities 

in  a  system  with  small  constant  shear  and  constant 

resistivity.   Wei-^  has  obtained  sufficient  conditions  for 

stability  in  the  special  case  of  constant  resistivity  and 

viscosity.   In  this  paper  we  investigate  the  exponential 

stability  of  static  equilibria.   Necessary  and  sufficient 

conditions  for  stability  are  obtained  without  any  of  the 

above  mentioned  approximations  or  limitations. 

We  adopt  the  model  of  reference  1,  with  the 
addition  of  an  arbitrary  nonconstant  viscosity.   The 
unperturbed  state  is  assumed  to  be  a  static  equilibrium 
solution  of  the  model  equations.   The  unperturbed  magnetic 
field  B^    (horizontal),  resistivity  r\^,    viscosity  v^,  and 
mass  density  p„  are  permitted  to  be  arbitrary  functions  of 
the  vertical  coordinate  z  consistent  with  the  equilibrium 
equations . 

Stability  is  discussed  in  terms  of  the  Lp-norms 
of  the  system  variables.   The  variable  ^(xft)  is  said  to 
be  exponentially  stable  if  and  only  if  for  every  e  >  0 
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there  exists  a  constant  M^  such  that  [J  |£,p  dv]  '^  =    ||4ll 

et 
<  M  e   for  all  t  _>  0.   Necessary  and  sufficient  conditions 

for  exponential  stability  of  the  sheet  pinch,  as  well  as 

the  maximal  growth  rate  of  the  unstable  pinch,  are  obtained 

directly  from  the  theorems  of  reference  8  and  Section  k   of 

this  paper.   These  theorems  extend  the  well-known  energy 

principle  for  conservative  systems  to  dissipative  systems 

of  a  rather  general  class,  and  are  valid  for  systems  with 

discrete  or  continuous  spectra.   The  maximal  growth  rate 

of  an  lins table  system  has  been  shown  to  be  the  supremum  of 

a  certain  functional  over  the  class  of  negative  potential 

Q 

energy  states.    Calculations  of  maximal  growth  rates  for 
the  unstable  pinch  will  be  presented  in  a  subsequent  paper. 
In  Section  2  we  state  the  nonlinear  model  equations, 
derive  the  equilibria  and  the  perturbation  equations,  and 
discuss  boundary  conditions.   Section  5  begins  with 
Theorem  5.1,  which  gives  the  necessary  and  sufficient 
conditions  for  the  exponential  stability  of  the  sheet 
pinch,  and  we  then  proceed  to  deduce  the  consequences. 
Section  4  provides  the  basis  for  Theorem  5.1.   Section  5 
leads,  in  particular,  to  the  following  conclusions  (the 
equilibrium  quantities  are  scaled  so  that  ^^(z)  is  fixed 
while  T],  the  magnitude  of  the  resistivity,  is  varied)! 
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1, 


The  stability  criteria  (but  not  growth  rates)  are 


independent  of  v„(z)  and  rj . 

dpQ 
2.   If  g  ^^  IS  positive  anywhere  (g  is  the  gravitational 

acceleration),  the  pinch  is  exponentially  unstable. 

This  holds  for  perfectly  conducting  or  insulating 

boundaries . 

5.   If  Bq  or  7]      is  constant,  stability  is  completely 

equivalent  to  gravitational  stability,  i.e.,  the 

pinch  is  exponentially  stable  if  and  only  if  g  — -  <   0 

dz  — 

This  holds  for  perfectly  conducting  or  insulating 
boundaries . 

4.  For  nonconstant  rj^,  the  system  is  always  unstable 
for  sufficiently  small  wavenumbers  with  insulating 
boundaries,  and  cannot  be  gravitationally  stabilized, 
but  it  can  be  gravitationally  stabilized  with  perfectly 
conducting  boundaries.   Strictly  speaking,  however, 
perfectly  conducting  boundaries  imply  zero  shear  in 
the  equilibrium  magnetic  field  —  see  Section  2(D). 

5.  Magnetic  stabilization  is  possible  with  perfectly 
conducting  boundary  conditions.   For  a  certain  class 
of  nonconstant  t^   (e.g.,  monotone)  and  no  gravity, 

the  imposition  of  a  properly  chosen  constant  zero-order 
magnetic  field  will  stabilize  a  previously  unstable 
system. 
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2.   The  Perturbation  Equations 

A.   The  Model 

We  consider  an  infinite  horizontal  layer  of  a 
magnetohydrodynamic  fluid  satisfying  the  following 
system  of  equations! 

V  •  V^  =  0  (2.1) 

|£  +  V''-  Vp  =  0  (2.2) 

p{||%  (v"- v)v^}  =  -  vp  -  pg e^  +  rx  r 


(2.3) 


+V'^(VV^)-  V^V^V-(VXV^)  XV  V 


•=•>     rr>.  ,  :::r>        -=•> 


E   +  V  X  B   =  Ti  J  (2.4) 

||"=  -  V  xr  (2.5) 

V  X  B^=  M-qJ'"  (2.6) 

V  .B^  =  0  (2.7) 

^  +  V":  vj,  =  0  (2.8) 
dt 


(2.9) 


The  equations  are  written  in  MKS  units*  M-q  is  the 
permeability  of  free  space.   The  quantity  p(x  ,t) 
denotes  the  mass  density,  V  (x  ,t)  the  fluid  velocity, 
B  (x  ,t)  the  magnetic  field,  P(x  ,t)  the  scalar  pressure. 
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J  (x  ,t)  the  electric  current  density,  v(x^,t)  the 
viscosity,  ri(x^,t)  the  resistivity,  g  the  gravitational 
acceleration,  and  E  (x  ,t)  the  electric  field.   We  choose 
a  right-handed  Cartesian  coordinate  system  (x,y,z);  z  is 


— > 


the  vertical  coordinate  with  corresponding  unit  vector  e"' 
The  gravitational  force  is  -pge"  .   Provided  the  r.h.s. 
of  Eq.  (2.9)  vanishes  for  the  zero-order  solutions  of 
Section  B,  its  precise  form  is  of  no  importance  in  the 
sequel  (see  Sections  B  and  C). 

The  quantities  E  and  J  may  be  eliminated  in 
favor  of  B^  and  V^  by  means  of  Eqs .  (2.4)  and  (2.6), 
resulting  in  the  set  of  equations  consisting  of 
Eqs.  (2.1),  (2.2),  (2.7),  (2.8),  (2.9)  and 

p{||  +  (V^V)V^]  =  -  VP  -  pge^  +  n-^(Vxt)x  B^ 

+  V^(vV^)-  V^'V^V  -(Vxv^)xVv 


(2. 3') 


||  =  -  V  X  (tiV  X  B^)PL-^  +  VX  (V^  X  B^)  .  (2.10) 

E  and  J  are  obtained  from  V  and  B  by  means  of 

J^=  tiQ-'-V  X  B^  (2.6') 

E^=  M'Q-'-TiV  X  b"^  -  V^x  B^  (2.4') 

With  the  exception  of  our  inclusion  of  a  nonconstant 
viscosity  v,  our  model  is  identical  to  that  considered 
in  Ref .  1;  the  reader  desiring  further  commentary  on 
these  equations  is  referred  thereto. 
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B.   The  Equilibrium 

The  houndaries  of  the  fluid  are  taken  to  be  z  =  0 
and  z  =  a.   The  equilibrium  quantities  (distinguished  by 
a  0-subscript)  are  assumed  to  be  independent  of  t  and  the 
horizontal  coordinates  x  and  yj  we  take 

Vq  =  0  (2.11) 

(no  mass  flow)  and  assume 

With  these  limitations,  the  equilibriiom  quantities  are 
allowed  to  be  arbitrary  functions  of  the  vertical  coordinate 
z  consistent  with  the  equations  of  Section  A.   Thus 
p^(z),  ^(-.(z),  and  V„(z)  are  permitted  to  be  arbitrarily 
prescribed  functions  of  z  on  [0,a]  such  that  p„  >  0, 
T]^  >   0,    and  Vq  >   0,    and  it  is  easily  seen  that  the  equations 
of  Section  A  will  then  all  be  satisfied  if  and  only  if 
the  remaining  equilibrium  quantities  £„,  T„j  Bq,  and  P„ 
satisfy  the  following  equations! 


^0  =  ^Ox^x  +  ^Qy^y  =  constant  (2.15) 

(E„  and  E„  are  arbitrarily  prescribed  constants) 

jj(z)  =  t1o^(z)eJ  (2.14) 

z 
bJ(z)  =  bJ(0)  +  ^XqEJ  X  e^  J  Ti-l(u)  du  (2.15) 

0 
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where  Bp^(O)  is  an  arbitrarily  prescribed  (constant) 

horizontal  field,  and  P^  is  given  by 

z 
Pq(z)  =  Pq(0)  +  J  [(Jq  XBQ).e^  -  spq]  du         (2.16) 

0 

where  Pr)(0)  is  an  arbitrary  constant.   This  completes 

the  description  of  our  equilibriiom  state,  which  is 

referred  to  as  the  "standard  case"  in  Ref.  1. 


G.   The  Linearized  Equations 

The  linearization  of  Eqs.  (2.1),  (2.2),  (2.^'), 
(2.7),  (2.8),  (2.9),  and  (2.10)  about  the  equilibrium 
solution  of  Section  B  leads  respectively  to  the  following 
equations  for  the  perturbed  quantities  (denoted  by  subscript 

1): 

V.  vj  =  0  (2.17) 


^  +  V,^p;(z)  =  0     (..^)  (2.18) 

hv^  — >  _  _i     _    _ 

Po  ^  =  -  ^^1  -  Pi<  +  V(^^^i)  ^^S  +  V^^^^o)^  ^1 

(2.19) 


z 


at 
and 
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V.  Bj^  =  0  (2.20) 

i  =  ?  (2.22) 


_i  =  _  VX(TiQVXBptx--'  -  Hq^Vx  (ti^VxBq)  +  Vx  (V^X^^.  (2.25) 
u  u 

Note  that  v   does  not  appear  in  Eqs .  (2^17 )- ( 2. 21 )  since 

Vj  =  0. 

The  perturbed  quantities  are  Fourier  analyzed  in 

the  horizontal  plane,  i.e.,  we  write  each  perturbed  quantity 

<t>-,  (x  ,t)  in  the  form 

i(k  X  +  k  v) 
<t.^(x^,t)  =  <f(z,t)  e   ^    ^       .  (2.24) 

Note  the  change  in  notation!   A  variable  (p,V  ,ri,P,J  _,B  ) 
without  a  numerical  subscript  shall  henceforth  denote  the 
transform  (a  function  of  z  and  t  only)  of  the  associated 
perturbed  variable  ( p,  ,¥-.  ,t]   p  T  ,  B-,   respectively). 
The  z-component  of  the  curl  of  Eq.  (2.19)  gives 

PqR  +  L-j_R  -  1[IqF   J^  =  q  B^  (2.25) 

Where  R  .  ||  ,  R-  .  ||  , 

_    -i(k  x+k  v) 
R(z,t)  ^  (V  XV^)^e    X   r   ^  ik^Vy(z,t)  -  ikyV^(z,t)  (2.26) 

J,(z,t)  =  ^.-^(V  XBp^  e  '   ^"""V   ^  ^-l(ik^By-  ikyB^)  (2.27) 

F(z)   =   ^x-l[k^BQ^(z)  +k^BQ^(z)]   =   ^^o^k^-Bj  (2.28) 

q(z)   =   i  Ho^(k^BQy  -  kyB^^)  (2.29) 

L^     =      -   DVqD  +  k^VQ  ,    D   =   ^  ,    k^  ^  kj  +  k^  .  (2.30) 

We  define  the  z-displacement  ^(z,t)  by 


t 

^(z,t)  =      J  V^(z,u)  du  +  ^qIz)  (2.51) 

0 

where  ^q(z)  is  an  arbitrary  fixnction  of  z.   The  integration 

of  Eqs.  (2.18)  and  (2.21)  with  respect  to  t  leads  to 

p(z,t)  =  p(z,0)  +  PqCq  -  Pq^  (2.52) 

Ti(z,t)  =  ri(z,0)  +  Tq^q  -  tIq^  (2.33) 

The  z-component  of  the  curl  of  Eq.  (2.25)  gives 


where 

L2  =  -  DtIqD  +  k^T^Q  .  (2.35) 

The  definition  of  R  and  Eq.  (2.I7)  yields 

V^(z,t)  =  k-2{ik^^'  +ikyR]  (2.36) 

Vy(z,t)  =  k'^{iky^'  -  ik^R]  (2.37) 


and  Eqs.  (2.27)  and  (2.20)  give 

B^(z,t)  =  k-2{ik^B^  +ikyHQj^)  (2.38) 

By(z,t)  =  k-2{ikyB^  -i^xVz^  (^-^9) 

Given  the  two  variables  C(z,t)  and  B  (z,t),  J  and  R 
are  determined  by  Eqs.  (2.25)  and  (2.34),  so  that 
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(2.40) 


V  ,  V  ,  B  and  B  can  then  be  obtained  from  Eqs .  (2.36)- 

(2.39). 

We  now  proceed  to  obtain  the  equations  determining 
the  quantities  ^(z,t)  and  B  (z,t).   We  multiply  the 
x-Gomponent  of  Eq.  (2.19)  by  ik  and  the  y-component  of 
Eq.  (2.19)  by  ik  ,  add  the  resulting  equations,  and  by 
means  of  Eqs.  (2.I7)  and  (2.20)  obtain 

k^P  =  -  Pq^  +  (Vq^  )  -  k^v^e  -  Vq^  +  k^VQ^  -  v^i 

-   IF'B^  +iFB;  -  k^^-V.Bj 

Substituting  this  expression  for  P  into  the  z-component 
of  Eq.  (2.19)  and  using  Eq.  (2.^2)  we  find 

L^i   +  L^k   -    gk^pQ^  +  iF[B^-  k\]  -  iF"B^ 

=  -  gk^[p(z,0)  +  Pq^q] 

where 

L^  H  -  DPqD  +  k^pQ  (2.42) 

L^  =  D^VqD^  -  2k^DVQD  +  k^VQ  +  k^VQ  (2.43) 

The  z-component  of  Eq.  (2.25)  and  Eqs.  (2.I7)  and  (2.33) 
imply 

(2.44) 

=    -  i  F'T1-1[T1(Z,0)   +  T)q^q] 

where 

-  10  - 


(2.41) 


-1    ?     ? 


(2.^5: 


The  term  iF(B^  "^^^z ^  "  ~^^^0^5^z  ^"  eliminated  from 
Eq.  {2. Hi)    by  means  of  Eq ,  (2.4-4),  and  we  write  the 
resulting  equation  together  with  Eq.  (2.44)  as  the 


single  matrix  equation 


(  L    0 


0    0 


2  2 


.  B 


^0 


^0^^  N  /  ^  X 


0 


'0 

-1 

'0 


\     z/ 


+ 


/  ^^o^F'Ti 


+ 


0 


''0 
iF'Ti 


-  gk  p 


2  I 
0 


0 


-  iF" 


L, 


'0 


/ 


(2.46) 


/-gk2[p(z,0)+pU„]+lx  Ti-lpF'  [Ti(z,0)+ii:?„] 


0^0- 


-iF'ii-l[Ti(z,0)  +  Til^„] 


'0 


'0"0- 


Eqs.  (2.15)   (2.15),  and  (2.28)  imply 

F'T-  =  ^x"-'-^k  b'  ii_  +  k  B'  T]^]  =  k  E„   -  k  E„  =  const.  (2.47) 
'0    O'-xQx'O    yOy'D-"     xOy    yOx  ^    '  ' 


so  that 


0  =  (F'Ti  )•  =  F"ii_  +  F'Til   <^=>  F"  ==  - 


'0' 


'0 


0 


„  I   I 

^0 


(2.48) 


This  completes  the  derivation  of  the  equations  for 
the  perturbed  quantities.   The  solution  of  Eqs.  (2.I7)- 
(2.23)  is  accomplished  by  solving  Eq.  (2.46)  for  ^  and  B^', 
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the  remaining  variables  are  then  obtained  from  Eqs .  (2.25), 
(2.32)-(2.34),  and  Eqs.  (2.^6 ) - (2.4o) . 

D.   Boundary  Conditions 

The  boundary  conditions  commonly  employed  in  the 

literature  are  those  associated  with  rigid  perfectly 

conducting  walls.   These  lead  to  the  simplest  boundary 

conditions  on  B  ,    namely  that  B  (O, t)  =  B  (a^t)  =  0. 

However,  perfect  conductors  require  that  E,    =0  thereon, 
^  ^  tan  ^ 

so  that  Eqs.  (2.I3)  and  (2.I5)  would  require  E^  s  0  and 
B(-.(z)  =  Bj^(O)  =  constant,  i.e.,  the  assumption  of  perfectly 
conducting  boundaries  implies  zero  shear  in  the  equilibrium 

magnetic  field  Bq(z).   To  avoid  this  specialization  we  shall 

■it- 
assume  rigid  perfectly  insulating  walls  at  z  =  0  and  z  =  a. 

Thus  we  require  that 

J^(0,t)  =  J^(a,t)  =  0  (2.49) 

V^(0,t)  =  V^(a,t)  =  0  (2.50) 

For  V„  >  0,  we  also  require  that 

V^(0,t)  =  Vy(0,t)  =  V^(a,t)  =  Vy(a,t)  =  0         (2.5I) 

which  together  with  Eq.  (2.I7)  implies 

V^(0,t)  =  V^(a,t)  =  0  (2.52) 


Perfectly  conducting  boundary  conditi^ons,  imposed  on 
the  perturbed  quantities  only  (i.e.,  Eq  ^   O),  will 
nevertheless  formally  be  considered  in  Section  5,  due 
to  their  common  usage  in  the  literature. 
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Eqs.  (2.26)  and  (2.5I)  lead  to 

R(0,t)  =  R(a,t)  =  0  (2.55) 

while  Eqs.  (2.31),  (2.5O),  and  (2.52)  give 

^(0,t)  =  |(a,t)  =  0  (2.5^) 

and 

|'(0,t)  =  l'(a,t)  =  0  (2.55) 

We  shall  require  that  the  arbitrary  function  ^(^{z) 
appearing  in  Eq.  (2.^1)  satisfy 

^0(0)  =  ?o(a)  ^  0  (2.56) 

so  that  Eqs.  (2.^1)  and  (2.50)  yield 

e(0,t)  =  ^(a,t)  =  0  .  (2.57) 

It  remains  to  obtain  the  boundary  conditions  on 
B  (z,t).   We  assume  that  Eqs.  (2.6)  and  (2.7)  hold 
everywhere  and  that  there  are  no  surface  currents  at 
the  walls.  These  conditions  then  imply  that  B  and  B. 
are  continuous  at  the  walls,  i.e.. 


B^(0+,t)  =  B^(0-,t),    By(0+,t)  =  By(0-,t) 


B^(a+,t)  =  B^(a-,t),    By(a+,t)  =  By(a-,t) 


(2.58) 


and 


B^(0+,t)  =  B^(0-,t),    B^(a+,t)  =  B^(a-,t)         (2.59) 


Eqs.  (2.20)  and  (2.58)  lead  to 

-  13  - 


B^(0+,t)    =    B^(0-,t),         B^(a+,t)    =   B^(a-,t)  (2.60) 

Thus  both  B     and   b'    are   continuous   at   z   =   0  and   z   =  a. 
z  z 

Outside   the  fluid,    the  equations  V  •  B-,    =   0  and 
V   X  B^   =   M-qJ^  =   0  imply   that 

->       ~^  ? 

Bj   =  V«J>^    ,        '^   *1  =0^         z<0,    z>a  (2.61) 


and  since  we  have  Fourier  transformed  in  the  horizontal 

^  i(k^x+k  y) 

plane,  we  set  <t>-.{x   ,t)   =   <ti(z,t)e       ^        and  obtain 

from  the  last  of  Eqs.  (2.6l) 


<l>"(z,t)  -  k^<l'(z,t)  =  0,      z<0,  z>a        (2.62) 
with  solution  <t>(z,t)  -  h(t)e*^.   Therefore 


-kh^(t)e 


-kz 


,   z  >  a 


B^(z,t)  =  J  (2.65) 

kh2(t)e^^   ,   z  <  0 

where  k  =  yk^ +  k^  =  v^  ,  so  that 

B^(a+,t)  =  -kB^(a+,t),   B^(0-,t)  =  kB2(0-,t)       (2.64) 

From  the  continuity  of  B  and  b'  at  0  and  a  we  obtain 
the  boundary  conditions 

B^(a-,t)  +  kB^(a-,t)  =  0,    B^ ( 0+, t )-kB^ ( 0+, t )  =0  (2.65) 

The  magnetic  field  B  (z,t)  outside  the  fluid  is  given  by 

B^(a,t)e-^(^-^),   z  >  a 

^z(^'^)  -  ^^  ,^,,  kz  ,       (2.66) 

B^(0,t)e       ,   z  ^  0  . 
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We  have  completed  the  derivation  of  the  necessary 
boundary  conditions,  and  are  now  in  a  position  to  specify 
the  domains  D^  of  the  operators  L^  (n  = 1,2,3,4, 5)  defined 
previously.   We  assume  henceforth  that  Pq(z),  t]Jz),    and 
Vq(z)   are  all  in  C  [0,a]  (the  class  of  all  twice 
continuously  dif f erentiable  functions  on  [0,a]).   In 
view  of  boundary  conditions  (2.53),  (2.49),  (2.54),  (2.55) 
and  (2.65)  we  take 

°1  "  °2  "  °3  ^  {f(z)|  f  e  c2[0,a],  f(o)=f(a)  =o}  (2.67) 
D^  ^  {f(z)  I  f  e  G^[0,a],  f  (  O)  =  f  (a  )  =  f  •  (  O)  =  f '  (a  )  =o]  (2.68) 
D^  ^  {f(z)  I  f  6  C^[0,a],  f'(a)  +kf(a)  =  0  =  f>(0)-  kf(0))(2.69) 

The  operators  L^  are  all  symmetric  on  D  and  in  particular 
L-j^  >;  0,  Lg  >  0,  L^  >  0,  L  ^  0  and  we  shall  assume  that 
L^^  ^  0  (it  suffices  for  example  that  v^  ^  0  on  [0,a]). 

E.   Dimensionless  Variables 

We  proceed  to  write  Eq.  (2.46)  in  dimensionless 
form.  To  this  end,  we  introduce  the  following  definitions: 

^  =  a'-'-z  ,       0  ^  z  <  a        (2.70) 

Ti   =     max   T    (z),         v   =     max   v    (z)    ,         p   =     max    pJz)  (2.7I) 

[0,a]    ^  [0,a]    "  [0,a]    " 

f(0    =  ^"^o(^^)'      h^'^)    "  r^PQ(aq,      s(C)   ^  v-^VQ(aC)      (2.72) 
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j'**^  T]   ^Eq      (constant)  (2-75) 

and   assuming  J*    =    JT  *  |    >   0, 

bM-qJ  "^  0 

r-^^.  ej  .  i^  =  -^  (2.75) 

(2.76) 

—  —  ^  t 

1c^  =  k%,      K   =  ka,        Y  =  ^^ — Q   J        u  =  7p 

W(C,u)      =      a"^^(a^,    uT)  (2.77) 

^(C,u)      H      (aHQj*)-l   B^(aC,uT)    .  (2.78) 

These  definitions  are  valid  for  J*  ^  0  (nonzero 
magnetic  shear)  which  represents  the  most  interesting 
case.   Indeed,  j  *  =  0  implies  that   B„  and  F  are  constants, 
and  it  will  be  shown  in  the  sequel  that  exponential 
stability  of  a  zero  shear  system  is  equivalent  to 
gravitational  stability,  i.e.,  the  system  is  exponentially 
stable  if  and  only  if  Pq  jl  0  on  [0,a].   Thus  the 
equilibrium  magnetic  field  E^  can  affect  exponential 
stability  only  if  it  possesses  nonzero  shear  (j*  ^^   O) . 

Equations  (2.l4),  (2.72),  and  (2.75)  imply 

Jq^^^  -   f'^f)j"*  (2.79) 
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Note  that  the  equilibrium  quantities  have  been 
scaled  so  that  by  keeping  J   fixed  while  letting  r|  vary, 
J^,  B^,  and  b  (C)  remain  fixed  independent  of  rj. 

Equation  (2. •46)  can  now  be  expressed  in  terms  of 
the  dimensionless  two-vector  0,  defined  by 

/W(C,u)\ 
0(C,u)  =  f    ^   J  (2.80) 

and  takes  the  form 

P©  +  K©  +  H©  =  G  (2.81) 

where  ©  =  ^r--  , 
ou 


(2.82) 


i-F".r 


(2.85) 


(2.84) 

(We  have  used  equation  (2.48)  to  put  H  in  symmetric  form.) 

r  =  (e^Xep.rf'f-2    (■  ^  ^)  (2.85) 

G  =  'c2[(e^-  b")r  -  Yh'(C)]  (2.86) 
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_     ..d^    d^    p^2  d    d     4     2  „, 


"dr     dc 


— > 

—>  K  f. 


V 


a  p 


aV(^tio)^/2 


(2.87) 
(2.88) 
(2.89) 


fGW(C,0)+K2{(e^r)(e^x^p-rT,(aC,0)f-2^-l-Yp(aC,0)ti^ 


\ 


,-2/— >w— >\  — >  I 


A-1 


V 


-iK:f-^(e;>€p-?^[f'W(C,0)+r"Ti(aC,0)] 


(2.90) 


/ 


The  domains  D  of  the  operators  L   (n  =  J>,^,5)    are. 


n 


n 


from  definitions  (2.67),  (2.68),  and  (2.69), 


D^  =  {f (C)  I  f  €  C^[0,1],  f(0)=  f(l)  =0} 


(2.91) 


D^  =  {f(C)|  f  e  C^[0,1],  f(0)=f(l)=f'(0)=f'(l)=0}  (2.92) 
D^  =  {f(C)  I  f  e  C^[0,l],f'(l)+*:f(l)  =  O=f'(0)-«:f(0))  (2.93) 

For  each  fixed  u,  0(^,u)  regarded  as  a  function  of 
C  is  an  element  of  the  Hilbert  space  /f  =   L^LO^l ] x  L^LO,! ], 
and  the  domains  D^,  D^,  and  D„  of  the  linear  operators  P, 

r         IS.  n 

K,  and  H  are  the  following  subsets  of  /r  '. 


Dp  =  D^  X  L2[0,l] 


(2.94) 


°K=' 


D,,    X  L„[0,1]    ,  vs    >    0  on    [0,1] 

^  ^  ^  (2.95) 


yf 


V    =   0 


Djj   =      L2[0,l]    X  D^ 


(2.96) 
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The  inner  product  ^ , ^  on  /r  is  of  course  given  by 

<Q-^,9^>   =    (W-L^W^)  +  (^1,^2)  (2.97) 

with  (W-L^W^)  =  J'^^1^2   ^^' 

and  it  is  easily  seen  that  P,  K,  and  H  are  synunetric 

on  their  domains  of  definition.   Furthermore,  we  have 

1 
<9,P©>=  JhClw'l^  +  «:^|w|^]  dC   >0  (2.98) 

0 
and 

1 

<©,Ke>  =  (w,L^w)  +  (n/ff)  J  f~^|(^^.b^)w  +  1^1  ^dC  (2.99) 

0 

so  that  Lji,  >  0  implies  that  K  ^  0. 
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3.   stability 

It  will  be  shown  in  Section  4  that  the  sheet 
pinch  is  exponentially  stable  to  disturbances  with  wave 
number  k"  =  T^a"-^  if  and  only  if  H  =  E{'ic^)    >_   0  for  that 
value  of  Ic^ ,    so  that  the  pinch  is  exponentially  stable 
to  arbitrary  disturbances  if  and  only  if  H  >^  0  for 
arbitrary  ''c'".   Necessary  and  sufficient  conditions 
that  H  be  positive  semidefinite  will  now  be  given,  and 
several  conclusions  regarding  stability  will  be  drawn 

A 

therefrom.   Since  H  is  independent  of  Vq  and  ri,  the 
stability  criteria  (but  not  the  growth  rates)  will 
of  course  be  independent  of  v„  and  t)  (-q  >  0). 

We  define  the  functional  F^[f]   for  all  real 

•v. 

^(C)    e   D       and   all    e  >   0  such  that   Q^  =   G+  e  >   0 
on    [0,1]   by 

1  2 

FgL^]   =  J   iirf  +  'c^[l-  ^]V'^3   d^  +  'c[^2(0)+/(l)]    .         (5.1) 

0  ^ 

Theorem  3.1.   HCk^)  >^  0  on  D„  if  and  only  if 
G  >  0  on  [0,1]  and  Fp[^]  >  0  for  all  real  Tp  e   Ti^   and 
all  e  >  0. 

Proof.   Our  previous  assumption  that  p„  e  C^[0,a] 
and  tIq  e  C^[0,a]  implies  that  G  =  G(k'',C)  e  C[0,1]. 
The  necessity  that  G  _>  0  on  [0,1]  is  obvious.   Suppose 
then  that  G  >  0  on  [0,1].   Let  W  =  W-^^  +  iW^  €  L2[0,l], 
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0 


'^  T  7 

Tp  =   Tp^   +  iTp^   e   Dp-;,  and  0  =  (  ,).   Then  for  e  >  0, 

1 
<9,He>  =  J  {W[GW+  licrf]    +  ^[-i/crW  +  L   f]}    dC 
0 

0 

0  ^ 

+  J  {g  JW^  +  -jj—  ]^  -  eW^  }  dC  +  F  J^-L  ] 
0  ^ 


where   R^[e,W,^]  =   f"""  {g^[W  +  '^   ]^   -    eW^}  dC  . 
±.  "^0       ^  ^e 

<\* 
Now  suppose  there  exists  a  real  V'olQ  ^  Dj-  such  that 

FgE^o-'  =  ^  ^  O3  fo^  some  e  >  0.   Then  there  exists  a 

real  W-,(C)  e  G   [0,1]  which  vanishes  identically  outside 

some  closed  subinterval  of  (0,1)  such  that 

I  G^[W^ G~^J  ^^   ^  "^  •   Choosing  V;^  =  0  and  W^  =  0, 

A 
equation  (3.2)  leads  to  <©,H©>  <  p  <  0.   Thus  H  >^  0 

implies  G  ^  0  on  [0,1]  and  F^[^]  _>  0  for  all  e  >  0 

and   all   real    ^  €    D^-. 

^  /W^+iW^^ 

Now    suppose  <e,H0>  =  A   <    0  for   some    ©  =[ 

^  \  ip  +lTp 

e   L2[0,l]   X   D    ,      Then  by   equation    (5.2),  ^ 
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and  this  holds  for  all  e  >  0.   Clearly  for  some 
sufficiently  small  e  >  0  we  have 

so  that  equation  (3.2)  yields 

Fe[^l]  +  F^li^^]    <  A  <  0  , 

which  implies  that  F  [^]  <  0  for  f  =   f     or   ip  =   ip^. 

Hence  G  >  0  on  [0,1]  and  F_[^]  >  0  for  all  real  ^  e  D^ 
—  fc.    —  ^ 

and  all  e  >  0  implies  that  H  >^  0.   This  completes  the 
proof. 

We  note  that  all  the  preceding  statements  of  this 
section  remain  valid  for  the  boundary  conditions 
^(0)  =  ^(1)  =  0  corresponding  to  perfectly  conducting 
boundaries,  provided  only  that  D^-  is  replaced  by 
D^  =  {V'l  ^ec^[0,l],  ^(0)  =  ^(1)  =  O}. 

The  quantity  r,  defined  in  equation  (2.85), 
vanishes  identically  in  C  if  f(C)  is  constant  or  if 
Ic      is  parallel  to  J^  .   In  this  case, 

»=( 

\   0       L^ 
and 

<e,H©>  =  -  /c2Y(W,h'W)  +  {4',L   f)    . 
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(3.3) 


Since  L,-  ^  0  on  D^-  and  D^-,    H  >^  0  if  and  only  if 
h'(C)  2  0  on  [0,1],  for  K^T  >  0.   Note  that  if  J*  =  0 
(i.e.,  Jq  =  O),  then  F'(z)  =   0  and  the  operator 
corresponding  to  H  in  equation  (2.46)  reduces  to 


2  ' 
0 


and,  as  above,  would  be  nonnegative  definite  if  and 
only  if  Pq(z)  <  0  on  [0,a],  for  kg  >  0.   Thus  if 
either  Jq  =  0  or  tj   is  constant,  then  exponential 
stability  is  equivalent  to  gravitational  stability, 
i.e.,  the  sheet  pinch  is  exponentially  stable  (for 
all  k  )  if  and  only  if  gpQ(z)  ;<  0  on  [0,a],  and  this 
holds  for  perfectly  conducting  as  well  as  insulating 
boundaries.   If  gpQ(z)  >  0  for  some  z  e  [0,a],  the 
pinch  will  be  exponentially  unstable  to  disturbances 
ith  k  parallel  to  J^  (regardless  of  whether  J„  =   0 


w 


or  T)   is  constant). 

The  preceding  paragraph  demonstrated  the  necessity 
that  g  P(-)(z)  jf^  0  on  [0,a]  for  the  sheet  pinch  to  be  stable 
(for  arbitrary  k  ).   In  general,  with  the  exception  of 
the  cases  T^^  =   O  (no  shear)  and  constant  r]^,    it  is  not 
sufficient.  To  see  this,  we  write 

G  =  /c2[Q^(epf'f-2  -  Yh']  (3.^) 

where 
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^C^iK")    "  [e^-  b"(C)](e^Xep.e^  ,   0  <  C  <  1       (5.5) 
Let  p(C)  =    f~  (r)  dr,  so  that  by  equation  {2.7h), 

b"{C)  =  r  +  p(C)eJxr  (5.6) 

Inserting  this  into  equation  (5'5)j  we  find 

and  for  C  >  0, 

Q^(ep  =  p"^{(e;-[pFjxr  +^])2_(e^  ^f]  (5.8) 


--> 


It  is  immediately  apparent  that  if  p  is  not 
parallel  to  e".  x  "e  ,  i.e.,  if  b   is  nonunidirectional, 
then  for  each  fixed  C  e  [0,1  ]j  Q,^{e    )    assumes  both 
positive  and  negative  values  as  e   varies.   Thus  even 
in  the  circumstance  that  gp„(C)  =  0,  p  •  e  .   =/  0   and 

A  A 

f'(C)   7/  0  for  some  C  e  [0,1]  guarantee  that  G  <  0  at 

<"  -> 

the  point  ^  for  some  range  of  values  of  e^  —  the  system 

would  then  be  unstable  to  disturbances  possessing  those 
values  of  e"  .   In  other  words,  nonconstant  ^   constant  p^ 
(or  no  gravity),  and  Bq(z)  nonxinidirectional  imply 
exponential  instability,  and  this  holds  for  perfectly 
conducting  or  insulating  boundaries.   Indeed,  if  B^(z) 
is  nonunidirectional,  the  system  can  be  made  lonstable 
for  any  gpQ(z)  by  choosing  an  ^  (z)  with  sufficiently 
large  gradients.   Suppose  now  that  p  is  parallel  to 
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e.  X  e   and  that  gp^  =  0.   Let  e.  =  e^,  so  that 
e .  X  e     =  e    ,    and  put  p  ^   ^e    .      Then 

J      Z      X  X 

G  =  <^{e^^-   eJ)^[P  +p{C)]ff-2  (5.9) 

so  that  a  necessary  condition  for  stability  when  gp'  =  0 
is 

f'(^)[P  +  P(^)]  >  0  on  [0,1]  (3.10) 

The  function  p  +  p(C)  is  strictly  increasing  on  [0,1]; 

thus  a  value  of  p  will  exist  such  that  equation  (5. 10) 

holds  if  and  only  if  for  every  pair  ^-1  j  ^p  ^  [Ojl]  ^or 

which  f'(C-j^)  <  0  and  f  {^^)    >  0  we  have  C-^  <  ^2'      ^^^^ 

a  function  f(C)  must  attain  its  maximum  (=l)  at  C  =  0 

or  C  =  1,  and  if  f(Cn)  =  min  f(C)  (>0),  then  f'(C)  <  0 

[0,1] 
on  [0,Cq]  and  f'(C)  >  0  on  [Cq,1],  and  equation  (5. 10) 

will  hold  for  p  =  -p(Cq).   If  f '  >  0  on  [0,1],  then  any 

positive  p  will  do,  while  if  f  <^  0  on  [0,1],  then  any 

t 
P  <  -p(l)  will  suffice.   We  conclude  that  when  gp„  =   0, 

the  system  will  be  exponentially  unstable  (G  <  0  for  all 

/c  >  0)  for  any  t^  (z)  such  that  Tq(z-]  )  >   0  and  t]' (z„)  <  0 

where  0  <_  z,  ^  Zp  j<  a,  and  this  holds  for  perfectly 

conducting  or  insulating  bovmdaries . 

The  function  Qr(e^^)  is  a  bounded  continuous 

function  of  ^  and  e^  for  any  given  positive  f  e  C  [0,1]. 

Ttierefore  given  eny  positive  ^r){z)    £  C  [0,a]  and  any 
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B^(0),  one  can  always  make  G  >  0  on  [0^1]  for  arbitrary 
k  7^  0  by  taking  gpo  sufficiently  negative  on  [0^1]  (see 
equation  (5.4)).   This  raises  the  possibility  of 
"gravitationally"  stabilizing  the  sheet  pinch.   Also, 
we  have  just  seen  that  for  certain  functions  f(C)  we 
can  make  G  _>  0  even  in  the  absence  of  gravity  by  choosing 
(3    (B^(0))  appropriately,  which  raises  the  possibility 
of  magnetic  stabilization.   However,  according  to 
Theorem  5.1,  G  ^  0  is  only  necessary  for  stability^  the 
functionals  F  must  also  be  nonnegative  if  the  system 
is  to  be  stable.   We  now  assume  G  _>  0  and  investigate 
the  functionals  F  .   We  assume  also  that  t]   is  nonconstant 
and  that  J   >  0,  thereby  excluding  the  trivial  cases  of 
Jq  =  0  and  constant  rj^,  which  have  been  dealt  with 
previously. 

Case  I.   Insulating  Boundaries 

Set  e  =   K^   >   0,    TpQ  =  1   +   'cC(l-C)  e  D  ,  and  let 

n(C).^  =    .     _.  "^ >0  (5.11) 

G+  'c      (e^-  b  )r  -  Yh'  +  1 

Note  that  n  is  independent  of  /c,  and  that  for  Ic^ 
nonparallel  to  e"^,  we  have    n(C)  dC  >  0.   Then 

1  1 

F  ni^Q]   =    ^^  J    {(1-2C)2  +  f^]    dC  +  2/c  -  J  ni^l   dC         (5.12) 
0  0 

so  that 
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lim  F  if    ]    =  -  ■  n  dC  <  0    (/c^x  e^  ^  O)        (^.13) 

K   _>0+  /C^    *-•         ^  T 


0 


J 


Thus  for  any  n,    however  small,  the  sheet  pinch  will  be 
exponentially  unstable  to  disturbances  with  sufficiently 
small  wavenumber  k,  provided  k  is  not  parallel  to  j!^, 
t]q  is  nonconstant,  and  |  Jq|  >  0.   Gravitational  or  magnetic 
stabilization  is  therefore  not  possible  for  such  equilibria 
with  insulating  boundaries,  unless  some  mechanism  is 
introduced  which  bounds  positive  k  away  from  zero  (e.g., 
finite  bo\indaries  in  the  x  and  y  directions).   It  is  clear 
from  equation  (3-1)  that  the  pinch  will  be  exponentially 

stable  to  disturbances  with  wavenximber  k  satisfying 

2                r                r^  1 

K      >   sup  I  -^ — ^^ J  . 

(^    (e^-  b  )r  -  7h' 


Case  II.   Perfectly  Conducting  Boundaries. 

Every  real   V(C)  e  D  satisfies  f  (V )^  dC 
1  -"  *^n 

^  T^      \      V'  dC,  so  that  equation  (3.I)  gives 

1                    2 
F^m    >  J  {  ir^  +  /c^  _  ^  ^  ^ ]  /  dC 

Therefore  the  system  will  be  exponentially  stable  to 
disturbances  with  wave  number  k  provided 


2 
s^P  f  ,_>  _>"" ]  <  TT^  +  /c^  (3.14) 

and  will  be  stable  for  arbitrary  disturbances  if 
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^'2  ^     ->  -> 

ic,X        (e^-b^)r  -Yh' 


}  <  TT' 


(3.15) 


Gravitational  stabilization  is  clearly  always  possible* 
given  any  ti„  and  Bq(0),  equation  (5.I5)  can  be  satisfied 
by  taking  gpi  sufficiently  large  and  negative.   Magnetic 
stabilization  is  possible  for  suitable  (e.g.j  monotone)  t]  * 
sufficient  conditions  for  exponential  stability  are  that 
SPn  —  ^  '^^  [Oj,a],  B^(z)  be  linidirectional^  equation  (3. 10) 
holds^  and 


S(C)  <  TT^  . 


0  <  C  1  1, 


(3.16) 


where 


s(C)  - 


for  0  <  C  <  1. 


,-2 


f'(C) 


P  +  P(C) 


lim  f' 

V^-^0 


f'(0 


P+p(C) 


P+p(C)  ^  0 


.  p(Cq)  =  -  p 


(3.17) 
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h.      stability  Theorems 

In  this  section  we  prove  the  assertion  made 
previously:   The  pinch  is  exponentially  stable  to 
disturbances  of  wavenumber  k  =  «:  a~   if  and  only  if 
hC'c^)  >_   0.   This  is  accomplished  by  a  straightforward 
application  of  several  theorems,  which  we  state  and  prove 
in  an  abstract  settingj  they  are  therefore  by  no  means 
restricted  in  application  to  the  specific  problem  of  the 
MHD  sheet  pinch. 

We  consider  functions  f(t)_,  defined  for  t  ^  0, 
whose  range  is  contained  in  some  inner-product  space 
with  norm  and  inner  product  denoted  by  ||  ||  and  (  ,  ). 
Such  a  function  f(t)  is  said  to  be  exponentially  stable 
if  for  every  e  >  0  there  exists  a  constant  M  <  oo  such 
that  i|f(t)||  <   MgC^*  for  t  >  0  (  if  f(t)  is  real  or 
complex-valued,  the  appropriate  norm  is  of  course  the 
modulus  of  f).   The  two  theorems  to  follow  provide  the 
sufficient  conditions  for  exponential  stability. 

Theorem  4.1.   Let  P,  H  and  iA  be  symmetric  operators 
defined  in  an  innner-product  space,  and  suppose  that 
H  >^  0  and  5  m   inf  i>^?^^  >  0.   Then  if  f(t)  is 
exponentially  stable,  so  is  every  solution  ^{t)    of 
the  equation 

i>X+   [n  + A]y.=  f{t)  t>0  (4.1) 
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Proof:  ^   (X.^X)  =  (PX.X)  +  0C,PX) 

=    (f-  [H+A]X,)^)    +   (X,f  -  [H  +A]-X) 
=  2  Re  {f,X)    -    2{X,hX)    <  2  Re  (f,;^). 

Integrating  from  0  to  t  and  using  the  exponential 

stability  of  f(t)  and  the  Schwarz  inequality ,  we 

obtain 

t 

5|IXlP  <  {X,PX)   1   ^o^P/^)  +  2M^  J  e"^IIX(T)ll  dx 

0 
so  that 

^  1/2 

IIX(t)ll  1   g(t)  -   [a  +  p  J  e^'^IIXII  dx]  (4.2) 

0 
where  a  =  ^''^  {'Xq>'P/(^)   and  p  s  25"-'-Mg. 

Thus  llXII/g  1  1,  and  since  g(t)  =  (p/2)g"^  ||):ile^^, 

t  t 

g(t)  -g(0)  =  (|)  J  e^llXllg-^(u)  du  <  (|)  J  e^  du 

0  0 

=  2|[e"-ll  . 

From  equation  (^.2)  we  obtain 

IIX(t)ll  <  g(0)  +  p(2e)-^e^^  <  [g(0)  +  p(2e)-^]e^^   . 

Theorem  4.2.   Let  P,  K  and  H  be  nonnegative  symmetric 
operators  defined  in  an  inner-product  space^  and  let 
C(t)  be  any  solution  of  the  equation 

•  •        • 

PC  +  kC  +  hC  =  r(t)  ,     t  >  0  .  (4.3) 
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(A)     If  ||r(t)||  <  M,  ||r(t)||  <  N  for  t  ^  0,  where  M 
and  N  are  constants,  and  if  5  =  inf  ^ /gel  >    0,  then 
there  exist  positive  constants  A,  B,  C,  D  and  E 
such  that 

||^(t)||  <  At  +  B  ,  t  >  0,  (i4.4) 

(C,PC)  ^  Ct^  +  Dt  +  E,        t  ^  0.  (4.5) 


(B)    If  r(t)  and  r(t)  are  exponentially  stable  and 

^f   nv,-p  (^j  fg  P4-  aK+  H]l)  ^  „  „^^  ^,T      n   ^u 
xf  inf  -^ — ^ — r-c — r-c — '-  >   0  for  all  a  >  0,  then 


TTTTT 

^  and  (C,PC)  are  exponentially  stable. 


•     • 


•     • 


Proof:  (A)  ^  [{k,H)  +   (C,hC)}  =  (pC+hC,C)+(C,pC  +hC) 

=  (r  -KC,C)  +  (C,r-KC)  =  2  Re(r,C)  -  2(C,KC) 
<  2  Re  {r/q    =  |^  2  Re  (r,^  -  2  Re  (r,C)  . 

Integrating  from  0  to  t  and  using  the  Schwa rz  inequality, 

we  obtain 

t 
(C,pC)    +   (C,hC)    <a  +  2||r||||Cll  +  2  J   ||r||    ||C  II  du  {U  .6) 

0 
where   a  =    (Cq,PCq)    +   (Cq,hCq)    -   2   Re    (rQ,CQ).    Hence 

t 
5||C|P   <  a   +  2MIICII  +  2N  J    ||C(u)||  du   , 

0 

so  that 

t 

(IICII   -   m-^f   <   a5-l   +  M^b-^   +  2NB-^  J   ||C  11  du    . 

0 
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t 

Set  p(t)  =   J  IICll  du,  A^  =   ab~^    +  M^5"^,  B^  =  2N6"^. 
Taking  the  square  root  of  both  sides  of  the  above 
inequality,  we  obtain 

1/2         n 

IIC(t)||  =  p(t)  <  (A^  +  B^p)     +  M6-^  (4.7) 

so  that 

-1/2 
^  {2B-l(A^  +  B^p)l/23  ^  p(A-L  +  B^p)      <  1  +  Mb-^A-l/^ 


Integrating  from  0  to  t  yields 

1/2     B.         _.   -,  /p      1/2 
(A-L+  B-^p)      1  (^)(1  +  M5-^A--'/^)t  +  A^    , 


and  equation  (4.4)  follows  at  once  from  equation  (4.7) 
From  equations  (4.4)  and  (4.6)  we  find 

t 

(C^PC)  <  a  +  2M(At  +  B)  +  2N  f  (Au  +B)  du 

0 

which  implies  equation  (4.5). 


(B)    Let  e  >  0,  and  set  C(t)  =  e^^C(t).   Then  ^(t) 
satisfies 


p|  +  Kgl  +  Hg^  =   f (t)  =      r  e    ^^    ,      t  >  0, 


where   K^  =   2eP  +K  >   0,    H^.  =    e^P  +  eK+  H   >   0,    and 


-et        ,,      et      -et 

-Et    ^    /  II*  11    ,  ^11.  ii\_-et 


inf    {^,tiA)/{i,^)    >   0.      Now    llfll   =    llrlle"^^    <   M,e^^   e'^^^M, 


and    llfll   =    llr  -er||e    ^^    <    (||r||  +£||r||)e~^^  1  N      +   eM     so 
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that  by  (A),  ll^(t)  ||  ^  At  +B  for  t  ^  0.   Therefore 
llC(t)||  =  ll^lle^"^  <  (At+B)e^*   for  t  >  0,  which  implies 
that  C  is  exponentially  stable.   Since  r,  r,  and  C  are 
exponentially  stable,  given  any  e  >  0  there  exist 
constants  M^.,  Ng  and  Mg  such  that  ||r||  ^  M  e^  , 


:^^  and  m\   <   M  -^^ 
from  equation  (4.6)  that 


rll  <   N^e  ^,    and  IICll  <  M^e  ^  for  t  >  0.   It  follows 


t 
{k,Fk)    <  a  +   2MgMge^^^  +  2Ng.Mg  J  e^^  du 

0 

<  [a  +  2MgMg  +  e'^NgMgle^^"^ 


which  proves  the  exponential  stability  of  (C,PC). 

Exponential  stability  of  the  pinch  for  H  >^  0 
follows  easily  from  Theorems  4.1  and  4.2.   The  crucial 
equations  are  (2.25),  (2.54),  and  (2.8l)  (or  its 
equivalent  (2.46)).   These  are  differential  equations 
in  t  for  the  quantities  R,  J  ,  W  (or  i) ,    and  f   (or  B  )', 
the  remaining  variables  are  all  obtained  algebraically 
from  this  fundamental  set  and  its  derivatives  with 
respect  to  t  and  z.   We  show  that  H  >_  0  implies  the 
exponential  stability  of  R,  J  ,  W  and  ii   in  the  Lp  norm. 
Suppose  then  that  H  >^  0,  where  H  is  given  by  equation 
(2.84).   It  is  a  simple  matter  to  verify  that  P,  K,  and 
H,  given  by  equations  ( 2. 82)- (2.84 )  with  domains 
(2. 94)- (2.96),  satisfy  the  hypothesis  of  Theorem  4.2  (B); 
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C,  given  by  equation  (2.90),  is  independent  of  t  so  that 
G  and  G  are  exponentially  stable;  therefore  Theorem  4.2  (B) 
implies  that  the  solution  ©(t)  of  equation  (2.81)  is 
exponentially  stable,  as  is  (O^PQ).   This  immediately 
establishes  the  exponential  stability  of  W,  f,^,];i, 
and  therefore  of  I,  B  ,  ^,    and  ?  .   From  the  fact  that 
1^  II?'  11  <  11^'  II  for  II4'  II  >  0,  we  infer  that  ?'  is 
exponentially  stable.   Equations  (2.25)  and  (2.54)  can 
be  written  as  the  matrix  system 

PX  +  (H+A);C  =  f 


with 


R  \      /  qB, 


J,A   ^  "V-q^^-  H^  Cq[no^-n(-.o)-TiQeQ]}. 


d 


The  hypothesis  of  Theorem  4.1  is  satisfied,  and  we 
conclude  that  R  and  J  are  exponentially  stable. 

The  case  k  =  0  requires  separate  consideration. 
In  this  circumstance,  the  perturbed  quantities  are 
functions  of  z  and  t  only,  and  equations  (2.I7)- (2.21) 
and  (2.23)  and  boundary  conditions  (2.50)-(2.5l)  lead 
to  (we  omit  the  subscript  1  on  the  perturbed  variables) 
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V     s    0,      J      =    0,      B      =   constant 
z  ^         z  z 

p(z,t)    =    p(z,0),  T(z,t)    =   il(z,0) 

z 
P(z,t)    =   -g  J  p(u,0)    du   -    t^Q^B^*    Bq  +  constant 

Po\   +  Ll\   =   ^o'-V^z^      \(0,t)   =   V^(a,t)    =    0 

(Vq   >   0) 
Po^y   +  Vy  =-^0^^Ox\'      V'^^    =  Vy(a,t)    =    0 


0 


M-^J      =    -   B    J  M-.J      =   B    , 

Ox  y^  0  y  x^ 


%\   +  ^2\   =    (^B^)'^         \(°'^)    =    \(^^t)    =    0 


^XqE^    +  L2By   =    (nBQy)    ,         By(0,t)    =   By(a,t)    =    0    . 


The  boundary  conditions  imposed  on  B  and  B  correspond 
to  perturbations  with  zero  net  electric  current,  i.e., 
they  imply  and  are  implied  by  perturbations  satisfying 

J  (z,t)  dz  =  0.   The  operators  L,  and  L^  are  nonnegative 
for  the  given  boundary  data,  and  exponential  stability 
follows  from  Theorem  4.1. 

Having  established  the  sufficiency  of  H  >_  0  for 
exponential  stability,  we  turn  to  the  demonstration  of 
its  necessity.   We  give  two  proofs.   The  first  is  quite 
general  in  application  but  assumes  the  existence  of  solu- 
tions to  equation  (2.8l)  satisfying  certain  initial  data 
--  an  assumption  obviously  justified  (indeed  required) 
on  physical  grounds  and  reasonable  mathematically,  as 
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the  initial  data  and  the  equilibrium  quantities  can  be 

oo 
chosen  to  be  G   functions.   The  second  proof  makes  no 

assumptions  --  there  we  show  that  if  H  is  not  positive 

semidefinite^  equation  (2.8l)  with  G  «  0  possesses  a 

normal-mode  solution  ©(z,t)  =  ©(z)e    with  co  >  0. 

Note  that  it  is  sufficient  to  demonstrate  the  existence 

of  an  exponentially  growing  solution  to  the  homogeneous 

equation  (2.81)^  since  given  any  i„{z) ,    we  may  choose 

r]{z,0)    and  p(z,0)  such  that  G  s  o  (see  equation  (2.46)). 

Theorem  4.5.   Let  P,  K,  and  H  be  symmetric  operators 
defined  in  an  inner-product  space,  P  >;  0,  K  >  0,  and 
-oo  <  5  ^  inf  H  <  0.   Suppose  that  for  all  sufficiently 
small  positive  e  the  equation  PC  +  kC  +  hC  =  0  admits 
of  a  solution  C(t;e)  for  t  >^  0  satisfying  the  initial 
data  C(0;e)  =  Ti,  C(o;e)  -  sr)  e  Np,  where  Np  is  the 
nullspace  of  P,  C  e  DpH  Dj^  O  D^  for  t  >  0,  and 
(t],Hti)  <  0.   Then  there  exists  positive  constants  a  and 
p  such  that  ||C(t;a)||  >  P  ^^^   for  t  >  0. 

Proof:   Let  e  >^  0,  and  define  H^  =  e^P  +  eK  +  H  and 
Kg  s  2ep  +  K.   Then  (Ti,HQr|)  =  (t],Hti)  <  0,  and  since 
(T^HgT])  is  a  continuous  function  of  £,  there  exists 
a  >  0  such  that  (ti,H^ii)  <  0  and   ^(t;a)  exists.   Let 
^(t)  =  C(t;a)e'°'^.   Then  ^(t)  satisfies  the  equation 
P^  +  K^t  +  H^^  =  0,  and  we  have  ^^  =  C(0;a)  =t], 
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Iq  =  ^(0;a)  -  aC(0;a)  =  C(0;a)  -  an  €  Np.   Now 
^I(Lp^)  +  (^,H^^))  =  -2(^,K^|)  <  0,  so  that 

=  (ri,H  Ti)  <  0,    and  we  have  lk(t)  ||^  >  6"-'-(ti,H  ri )  >  0 

for  t  >  0.   Therefore  IIC(t;a)ll  =  ll^lle"'^  >  [6-1(ti,H  n)]^^e"^. 

Theorem  4.4.   Suppose  v^  =  0,  k:  >  0,  and  that  f(C)  and  h(C) 

are  positive  continuously  diff erentiable  functions  on  [0,1]. 

/q  ngy 
Let   inf  yg  Q\      ^  0.   Then  there  exists  a  nontrivial 

solution  to  the  equation 


P©  +  KG  +  H0(C,t)  =  0  (4.8) 


(the  homogeneous  form  of  equation  (2.8l))  of  the  form 

/W(C)\  ^^ 
©(C,t)  -I  )e   where  O  is  a  positive  constant. 

Proof.   Let  Jf^    denote  the  Hilbert  space  Lp[0, 1], 
9/  =  L^CO^l]  X  L2[0,l],  and  let  D  denote  either  D  or  D  . 
The  hypotheses  on  h  and  /c  guarantee  the  existence  of  the 
positive  compact  Hermitian  integral  operators  K^  and  K^  on 
Jf^     possessing  the  following  properties! 


t\j  t\. 


K^L^  =  I   on  D^  ,     K^L^  =1   on  D^  (4.9) 

3  5         3       5  5         5  V   ^/ 

L  K  =  I  on  C[0,1],  L  K  =  I  on  C[0,1]  (^.lO) 

K^Of)  c  C[0,1],      K^()/)  c  C[0,1]  (4.11) 
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Define  k  =  K^^ ,    ^^   =  K^   .   Then  k^  and  k^  are  positive 

compact  Hermitian  operators  on  Tr.   It  will  be  convenient 

to  denote  the  elements  of  the  symmetric  matrix  operators 

K  and  H  by  K.  .  and  H.  .   (i,j  =1,2).      For  each  real  03  we 
J- J      -L  J 

define  the  operators 


0)^1  0 


'k^  [^K^l+H-Li  ]k^    k^  [coK^2+%2  ^^5 


^    V  0   I  /      '^    Vk^LcDK^i+H^ilk^   k^^K22k5 

The  hypothesis  on  f  and  h  insures  that  H-j^-^^  H-^^^^  ^21 
as  well  as  all  the  K. .  are  continuous  functions  of  C 
on  [Ojl].   Thus  C^  is  compact.   P^  and  G^  are  clearly 

Hermitian.   Let  F(a))  s   mf  .   F(a))  is  a 

real-valued  continuous  function  of  cd  on  (-00,00).   It 

is  not  difficult  to  show  that  F(cjo)  is  nonnegative  for 

all  sufficiently  large  od,    and  that  inf  y'  qV   ^  0 

Dy  <^>^^ 
implies  F(o)  <  0.   Therefore  there  exists  Q  >  0  such 


) 


that  F(0)  =  0,  and  the  definition  of  F(0)  yields 

sup   g  p  ev  =1.   It  follows  from  well-known  theorems 

oil  compact  Hermitian  operators  that  there  exists  a 

/^.(C)\ 
nontrivial  i  =  [  \  &    \L     such  that  P^^  =  CU4.   Let 

W(C)  =  k  L   and      V'(C)  =  kj-^^.   Then  we  have 

O^W  =  -  K^{[OKj_-L+  H^-^JW  +  [OK-L2+  H-i_2]^}  (^.12) 

^  -  -  K^{[OK2i+  H2;^]W  +  [OK^g]^}  (^.13) 
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It  follows  at  once  from  (4.11)  that  W  and  V  are  continuous 

'X. 

on  [0,1];  (-^1,10)  then  implies  that  W  e  D-,  and  ^  e  D 

y  5* 

Multiplying  both  sides  of  equations  (4.12)  and  (4.1^)  on 
the  left  with  L^  and  L  ,  respectively,  we  obtain 

{O^P  +ok  +  h}(M=o, 
and  the  proof  is  complete. 
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